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ABSTRACT 

We report the detection of the radio and infrared counterparts of the ring nebula 
around the WN3(h) star HD 211564 (WR152), located to the southwest of the Hn 
region Sh2 132. Using radio continuum data from the Canadian Galactic Plane Survey, 
we identified the radio counterparts of the two concentric rings, of about 9' and 16' 
in radius, related to the star. After applying a filling factor / = 0.05-0.12, electron 
densities and ionized masses are in the range 10-16 cm~ 3 and 450-700 M Q , respectively. 
The analysis of the Hi gas emission distribution allowed the identification of 5900 M Q 
of neutral atomic gas with velocities between -52 and -43 km s _1 probably linked 
to the nebula. The region of the nebula is almost free of molecular gas. Only four 
small clumps were detected, with a total molecular mass of 790 M©. About 310 M© 
are related to a small infrared shell-like source linked to the inner ring, which is also 
detected in the MSX band A. An IRAS YSO candidate is detected in coincidence with 
the shell-like IR source. 

We suggest that the optical nebula and its neutral counterparts originated from 
the stellar winds from the WR star and its massive progenitor, and are evolving in 
the envelope of a slowly expanding shell centered at (I, b) — (102?30, -0?50), of about 
31 pc in radius. The bubble's energy conversion efficiency is in agreement with recent 
numerical analysis and with observational results. 
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1 INTRODUCTION 

The strong mass flow and high UV photon flux of Wolf- 
Rayet stars originate the so called ring nebulae, which can 
be identified as ring-like optical features in the environs of 
the stars. Searches for optical ring nebulae around WR stars 
performed in the last decades (e.g. Chu 1981; Heckathorn et 
al. 1982; Miller & Chu 1993; Marston et al. 1994, 1995), 
resulted in the detection of more than 60 ring nebulae, 70% 
of which are linked to WN stars. High angular resolution 
radio continuum observations led to the detection of the 
radio counterparts to 25% of the ring nebulae, allowing the 
estimation of several physical parameters (e.g. Cappa et al. 
2002). 
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In addition, spectral line observations of the neutral 
atomic and molecular gas helped to determine the charac- 
teristics of stellar wind bubbles, since they revealed large 
amounts of neutral gas linked to the nebulae, and allowed 
to include neutral gas in the estimate of the bubble's ener- 
getics (Cappa 2006). Conclusions from these studies, which 
indicated that a small amount of the stellar wind energy 
released by the central star is converted into kinetic energy 
of the bubbles, are compatible with both observational and 
theoretical studies (e.g. Chu et al. 1983; Oey 1996; Freyer et 
al. 2003, 2006; Cooper et al. 2004). 

In the last years, available infrared images obtained 
from the MSX Galactic Plane Survey (Price et al. 2001) 
and GLIMPSE survey (IRAC images, Benjamin et al. 2003) 
allowed the identification of interstellar bubbles in the near- 
and mid-IR, showing the presence of photodissociation re- 
gions at the interface between the ionized and molecular gas 
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Figure 1. Top left panel: DSS2 R image of the ring nebula around HD 211564. The cross marks the position of the WR star. The color 
scale is in arbitrary units. Top right panel: CGPS radio continuum image at 1420 MHz. Contour lines are 7.0 and 7.4 K. Bottom left 
panel: CGPS radio continuum image at 408 MHz. Contour lines are 62, 65, 68, 71, 74, 100, 150, and 200 K. Bottom right panel: Overlay 
of the optical (in color scale) and 1420 MHz (contour lines) images. 



(Churchwell et al. 2006, Watson et al. 2008, Watson et al. 
2009). 

Among the ring nebulae detected optically, the one 
around the WN star HD 211564 is particularly interesting. 
The morphology of this optical ring nebula can be described 
by two concentric ring structures having about 17' and 36' 
in diameter, detected in [OIII], [SII], and Ha + [Nil] lines 
(Heckathorn et al. 1982). The nebulae are easily identified 
in the DSS2 R image displayed in the upper left panel of 
Fig. 1. The symmetry of the inner ring is suggestive of a 
quite homogeneous interstellar medium, except possibly in 
the south-southwest quadrant. Both structures are brighter 
in [OIII] lines than in Ha + [Nil], following the tendency 
for nebulae around early WN stars (Heckathorn et al. 1982) . 
The inner structure is sharp in appearance and almost com- 
plete, with the WR star displaced 3'.7 northeast from the 



geometrical centre. On the contrary, the outer shell is more 
diffuse and less well defined. 

This ring nebula is located at the southwest of the Hn 
region Sh2-132, whose brighter section (not shown here) 
is excited by the WR star HD 211853, at (l,b) = (102?78, 
-0?65), and a few O-type stars. Although both HD 211853 
and HD 211564 are linked to Sh2-132, they appear related 
to quite distinct structures. 

Previous radio observations of Sh2-132 at 1420 MHz 
(Felli & Churchwell 1972) show a relatively large area of 
faint emission extending towards the region where the ring 
nebula is present. The inner shell was also detected in the 
far infrared by Marston (1996), who found an infrared shell 
of 21' in diameter surrounding the WR star. 

HD 211564 (ee LSIII+55 30 = WR152) is a WN3(h) 
star (Smith et al. 1996) located at (/,&) = (102?23, -0?89) 
or (RA,Dec[J2000]) = (22 h 16 m 24 s , +55°37'37"). The star is 
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considered a probable member of the Cep OBI association 
(Lundstron & Stenholm 1984) located in the Perseus spiral 
arm. 

The terminal wind velocity of HD 211564 is in the range 
1800-2100 km s" 1 (Rochowicz & Niedzielski 1995; Hamann 
& Koesterke 1998; Niedzielski & Skorzynski 2002), while a 
mean value for the mass loss rate of WN-w stars is about 
5xlO~ 6 M yr" 1 (Crowther 2007). 

In this paper, we report on the radio counterpart of the 
ring nebula around the WN star HD 211564 first identified 
by Heckathorn et al. (1982) and analyze the distribution of 
the ionized and neutral gas using radio continuum and Hl- 
21cm line data from the Canadian Galactic Plane Survey 
(CGPS, Taylor et al. 2003), CO data from the Five Col- 
lege Radio Astrophysical Observatory (FCRAO, Heyer et 
al. 1998) , and infrared data from the IRAS and MSX satel- 
lites. 



2 DATA SETS 

The analysis of the ionized and neutral atomic gas in the 
environs of HD 211564 was performed using data from the 
CGPS obtained with the Synthesis Telescope of the Domin- 
ion Radio Astrophysical Observatory (DRAO), in Canada. 
This telescope performed interferometric observations of the 
21-cm Hi spectral line, and, simultaneously, continuum emis- 
sion in two bands centered at 1420 MHz and 408 MHz. 
Single-dish data were incorporated into the interferometric 
images ensuring accurate representation of all structures to 
the largest scales. 

To investigate the Hi emission distribution, we ex- 
tracted a data cube centered at (l,b,v) = (102? 50, -1?0, 
-53.4 km s _1 ) and analyzed a region of l?5xl?5 in size. 
The synthesized beam is Ifl9x0.'98, the rms noise is 3 K in 
brightness temperature (7b), and the velocity resolution and 
channel separation are 1.32 and 0.824 km s , respectively. 
The Hi images were convolved to a 2. '5 x2.'5 beam size to 
facilitate the identification of structures. The observed ve- 
locities cover the range -165 to +57 km s _1 . 

The radio continuum data have synthesized beams of 
3!4x2.'8 and l.'0x0.'82 at 408 and 1420 MHz, respectively. 
The measured rms image noises are 1.1 and 0.063 K at 408 
and 1420 MHz, respectively. 

The 12 CO (1-0) line data at 115 GHz were obtained from 
Brunt & Heyer (2009). The angular resolution is 45", the 
velocity resolution 1.0 km s _ , and the rms noise, 0.7 K 
(main beam brightness temperature). 

Infrared images at different wavelengths were used to 
investigate the dust distribution in the region of the nebula. 
High resolution IRAS images (HIRES) at 12, 25, 60, and 
100 fim, were taken from IPACq The angular resolution is 
in the range l'-2'. 

Images at 8.28, 12.13, 14.65, and 21.3 (im (bands A, C, 
D and E, respectively) taken by the MSX satellite were also 
obtained from IPAC (18'.' 4 in angular resolution). 



1 IPAC is funded by NASA as part of the IRAS extended mission 
under contract to Jet Propulsion Laboratory (JPL) and California 
Institute of Technology (Caltech). 



Additional radio and optical images of the region were 
retrieved from the Skyview web pagqj. 

Finally, to investigate the presence of stellar formation 
activity in the outskirts of the nebula, we used infrared point 
sources from the MSX, 2MASS, and IRAS catalogues. 



3 RESULTS 

3.1 The ionized gas 

The radio emission at 1420 and 408 MHz, along with an 
overlay of the optical and radio emissions at 1420 MHz, 
are displayed in Fig. 1. Both the inner and outer features 
identified at optical wavelengths have well detected radio 
counterparts, as observed at 1420 MHz (upper right panel). 
The lower right panel shows the striking correlation between 
the optical and 1420 MHz radio continuum. The inner ring 
is not complete. However, very low level radio continuum 
emission is present near (l,b) = (102? 07,-1? 00). The outer 
ring is incomplete both at optical and radio wavelengths. 
On the other hand, at 408 MHz only the inner structure can 
be easily identified, while emission linked to the outer op- 
tical feature is discernible only near (l,b) = (102?0, -0?67) 
(Fig. 1, bottom left panel). The inner and outer rings, as 
estimated from the optical and radio images, are about 19' 
and 32' in size, respectively. 

The nebula is also detected at 2.7 GHz (Fiirst et al. 
1990) and at 4.85 GHz (Condon et al. 1994). Faint emission 
at 8.35 and 14.85 GHz (Langston et al. 2000) most probably 
associated with the ring nebula can also be identified. 

In order to estimate flux densities of both the inner 
and outer structures at 1420 MHz and to derive the radio 
continuum spectral index, we analyzed first the nature of 
the small-diameter radio sources projected onto the optical 
nebula. The coordinates of such sources are listed in Ta- 
ble 1, along with the measured flux densities at 408 and 
1420 MHz, the spectral index derived from the CGPS data, 
the catalogued flux densities at other wavelengths, and their 
identification. 

Sources 1, 2, 4, 5, and 6 appear projected onto the neb- 
ula. The location of sources 7 and 8 suggests that they are 
probably unconnected to the features we are interested in. 
Sources 1, 2, 5, 6, and 8 are clearly non-thermal in nature, 
probably extragalactic sources. 

After subtracting the contribution of the compact 
sources, and taking into account different background emis- 
sions, we estimated the flux density of the inner and outer 
features around HD 211564 from the image at 1420 MHz. 
Our estimates are S„ — 0.6±0.1 Jy and S„ — 1.2±0.1 Jy for 
the inner and outer structures, respectively. 

The CGPS images were used to evaluate the spectral 
index distribution of the radio continuum emission around 
HD 211564. As a first step, we removed the point sources 
from the region of interest and convolved the 1420 and 408 
MHz images to the same circular resolution of 3^4. We then 
constructed TT-plots, in which the brightness temperature 
Tb at one frequency is plotted point-by-point against the 
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Figure 2. TT-plots corresponding to (a) the entire region 
within an angular radius of 24' from the center of the 
inner ring, (b) the northwestern outer ring and (c) the northern 
half of the inner ring. For all plots, a background of 6 K and 57 K 
has been subtracted from the convolved images at 1420 and 408 
MHz, respectively. 



Figure 3. Top panel: IRAS (HIRES) image at 60 /im of the ring 
nebula in color scale and contours. Contours are 25, 30, 35, 50, 
70, and 90 MJy ster -1 . Bottom panel: Overlay of the images at 
60 fim (in color scale) and at 1420 MHz (in contours). 



brightness temperature at the other frequency. The tem- 
perature spectral index j3, where Tb oc is directly re- 
lated to the slope of the straight line fitted by regression, 
= log(slope)/log(1420/408). The usual flux density spec- 
tral index a (S v oc u a ) is simply a — f3 — 2. 

The results of the TT-plot analysis are presented in 
Figure 2(a-c). In each case, the analysis was limited to re- 
gions brighter than 6 K at 1420 MHz and 57 K at 408 MHz. 
The first panel (Fig. [2^,) is a TT-plot of the entire region 
within an angular radius of 24' from the centre of the in- 
ner ring at (l,b) = (102? 20, -0?92). Although a significant 
amount of scatter is present on the diagram, a convincing 
fit is obtained yielding a — 0.24 ± 0.05. Since this value 
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Table 1. Compact radio sources towards the ring nebula around WR152. 



No. 


(l,b) 


Sl420 


S408 


a 


S365 


S610 


S1420 


S4850 


Identification 




(°) 


mjy 


mJy 




mJy 


mJy 


mJy 


mJy 








(1) 


(1) 




(2) 


(3) 


(4) 


(5) 




1 


102.355,-1.065 


165±5 


409±22 


-0.7±0.2 


677+78 


245 




52±8 


87GB[BWE91] 2216+5517 


2 


102.308,-0.961 


187±6 


630+20 


-1.0+0.1 


725+56 


325 




91±12 


87GB,221520.4+552034 


3 


102.206,-0.160 


10±1 










9.7±0.6 




NVSS J221720+552309 


4 


102.122,-0.932 


5.8±0.5 










6.6±0.6 




NVSSJ 221556+553135 


5 


102.161,-0.896 


17±1 


35±4 


-0.6±0.3 






13.6±0.6 




NVSS J221601+553446 


6 


101.986,-0.856 


44±1 


87±4 


-0.6±0.1 










7C2213+5516 


7 


102.519,-0.743 


32±1 








40 


23.0±0.9 




NVSS J221731+555424 


8 


102.466,-0.619 


30±1 


66±7 


-0.6±0.3 




65 


19.2±0.8 




NVSS J221632+560118 



References: (1) this paper, derived using CGPS data; (2) Douglas (1996), Gregory (1991); (3) Harten et al. (1978); (4) Condon et al. 

(1998); (5) Becker (1991) 



is somewhat different from the canonical value of -0.1 for 
optically thin bremstrahhmg, we applied the same proce- 
dure to a number of sub-regions around HD 211564 to look 
for possible spatial variations. Figures [2jb) and (c) show 
the TT-plots corresponding to the northwestern outer ring 
near (l,b) = (102?08, -0?67) and the northern half of the in- 
ner ring, respectively. There seems indeed to be significant 
variations. As a further test, we evaluated the spectral in- 
dex of the nebulosities around nearby HD 211853 and found 
q = 0.00 ± 0.03 (Vasquez et al. 2009), a value clearly more 
consistent with optically thin bremstrahlung. 

3.2 Infrared emission 

The distribution of the infrared emission at 60 /im and its 
correlation with the radio emission at 1420 MHz is displayed 
in Fig. 3. The brighter sections of the inner shell are well de- 
tected in the far infrared. Emission from the outer feature 
is hardly detected. The far IR fluxes estimated for the ring 
nebula at 100 and 60 ^im are Fioo = 660+130 Jy and Feo 
— 340+100 Jy, respectively. The quoted uncertainty was de- 
rived by using different values for the background emission. 

Particularly interesting is the strong and extended IR 
source located at (l,b) — (102?20, -0?73). The source is de- 
tected in the four IRAS bands. Its ring morphology is shown 
in the upper panel of Fig. 4, which displays an overlay of the 
emission at 25 (color scale and white contours) and in 
the radio continuum at 1420 MHz (black contours). This IR 
source partially coincides with a strong radio emitting re- 
gion. The extended IR source has a clear counterpart in the 
MSX band A, although the emission is faint (Fig. 4, bottom 
panel). Note that the brightest region at 8.3 ^m does not 
emit at 1420 MHz. Measured IR fluxes at 100, 60, 25, and 12 
fim are 127+23, 53+3, 5.4+1.2, and 3.8+0.2 Jy, respectively. 
The estimated dust mass is 0.5+0.3 M©. 

Following the procedure described by Cichowolski et al. 
(2001), we derived the color temperature of the dust associ- 
ated with the ring nebula and the strong IR source. Taking 
into account different values for the background emission, 
we found similar dust color temperatures Td = 30+4 K for 
the ring nebula and the strong source. The range of tem- 
peratures corresponds to n = 1-2, where the parameter n is 
related to the dust absorption efficiency (k„ oc v n ). The 
dust temperature derived for both features is typical for Hn 
regions. The dust mass linked to the nebula is 2+1 Mq. 
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Figure 4. Top panel: Overlay of the emission at 25 (im (in color 
scale and white contours) and the emission at 1420 MHz (in black 
contours). The white contours are 10, 12, and 15 MJy ster -1 , 
while the black contours are 7.4, 8.0, and 8.5 K. Bottom panel: 
Overlay of the emission at 8.3 fim (MSX band A) (color scale) and 
the CO(1-0) emission (in contours) within the velocity interval 
from -49.4 to -46.4 km s — 1 (see Sect. 3.3). The contours are 3.0, 
7.0, 11.0, and 15.0 K km s" 1 . 



3.3 Neutral gas distribution 

The analysis of the Hi datacube indicated that gas probably 
linked to the ring nebula is present in the velocity interval 
from -57 to -37 km s _1 . Figure 5 displays overlays of the Hi 
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Table 2. IRAS point sources classified as YSO candidates. 













Fluxes 




No. 


(l,b) 


Name 


12^im 


25/im 


60/j.m 


100/im 






O 




Jy 


Jy 


Jy 


Jy 


L 


1 


102.122,-0.790 


22135+5523 


0.26 


0.30 


3.55 


22.20 


166 


2 


102.203,-0.723 


22137+5529 


0.91 


0.80 


20.75 


74.75 


705 


3 


102.440,-0.895 


22158+5529 


0.25 


0.30 


3.60 


21.60 


164 



emission within the velocity interval from -55.5 to -39.0 km 
s _1 (in color scale) and the radio continuum at 1420 MHz 
(contours). Each Hi image is the result of averaging the Hi 
emission within 3.3 km s _1 . 

The gallery of six Hi images was produced by removing 
the background in each channel map using DRAO software. 
The removed background was estimated as the average over 
the displayed image excluding a circular area of 28' in radius 
centered on the WR star. An inspection of the images shows 
that the gas distribution is quite complex, showing several 
arc-like features. 

The Hi emission in the range -47 to -39 km s _1 (Fig. 
5d, e, and f ) is characterized by the presence of a large cavity 
centered at (l,b) — (102?30, -0?50), encircled by an almost 
circular envelope. The structure is about 30' in radius. The 
inner and outer rings linked to the WR star appear projected 
onto the weakest section of this large shell. 

The images corresponding to the above mentioned ve- 
locity interval reveal the existence of neutral gas linked to 
the ring nebula. The presence of Hi emission encircling the 
inner shell is easily identified between -47 and -42 km s _1 
near (I, b) = (102? 35, -0?90) (Fig. 5d and e). This neutral gas 
forms a partial ring. The outer ring appears also bordered 
by Hi gas near (l,b) = (102?00, -0?60). Note the presence of 
gas projected inside the inner ring. 

Hi emission in the velocity range from -52 to -49 km 
s _1 (Fig. 5b and c) is seen projected all over the area cov- 
ered by the rings, being even observed inside the inner ring 
and in the region between the ionized rings. The partial Hi 
ring detected near (l,b) = (102?35, -0?90) (Fig. 5d and e) 
is also detected in this velocity interval. The emission ex- 
tends towards (l,b) = (102?08, -1?17), beyond the ionized 
arcs. Enhanced Hi emission bordering the inner ring is also 
present at -55.5 km s _1 (Fig. 5a). 

The complex Hi gas distribution between -47 and -39 
km s _1 suggests the existence of two structures with similar 
velocity ranges partially superpossed in the line of sight: the 
large Hi shell centered at (l,b) = (102? 30, -0?50) and the 
partial Hi rings linked to the ring nebula. 

The large Hi shell and its correlation with the optical 
emission is shown in Fig. 6. The upper panel displays an 
overlay of the radio continuum emission at 1420 MHz (in 
contours) and the Hi emission between -47 and -37 km s _1 
(in color scale). No background emission was removed 
from the Hi image displayed in this fugure. The mid- 
dle panel displays the DSS2-R image of the same area, and 
the bottom panel shows an overlay of the optical image and 
the Hi image of the upper panel (in contours). The satu- 
rated region at (Z,6) = (102? 80, -0?70) in the optical im- 
age corresponds to the main body of the Hn region Sh2-132 
(Vasquez et al. 2009) . The optical filament detected between 



(/,&) = (102? 40, -0?75) and (l,b) = (102? 17, -0?35) and its 
faint extension towards 6 > -0?35 defines an almost com- 
plete ring-like structure coincident with the inner border of 
the large shell (see the bottom panel). Note that the bright- 
est section of this filament corresponds to the outer ring. 
No background was removed from the Hi image displayed in 
this figure. 

To summarise, gas clearly linked to the ring nebula 
is detected between -52 and -43 km s~ x (Figs. 5b to 5e), 
while small patches of Hi gas are detected beyond these val- 
ues. The large Hi shell centered at (l,b) = (102?30, -0?50), 
of about 30' in radius, has an optical counterpart and is 
detected from -47 to -37 km s _1 in Hi emission. Neutral 
atomic gas belonging to these structures is partially coinci- 
dent in the line of sight. 

The distribution of the molecular gas was investigated 
using the 12 CO(1-0) datacube, in the velocity interval from 
-105 to +24 km s _1 . Within the area of interest, only a 
few patches of CO emission were detected within the ve- 
locity ranges from -63 to -47 km s~\ -30 to -28 km s~\ 
and -9 to -2 km s _1 . Molecular gas having velocities v < 
-47 km s _1 is probably linked to the ring nebula. The inte- 
grated emission between -56.5 and -46.4 km s~ is displayed 
in Fig. 7 superimposed onto the radio continuum emission. 
Only four molecular clumps are detected, and three of them 
coincide with the optical rings (clump A at [l,b] = [102? 13, 
-0?80], clump B at [l,b] = [102? 23, -0?70], and clump C at 
[l,b] = [102? 32, -0?96]). Note also that CO clump D at (l,b) 
= (102? 38, -0?85) and Clump B are projected onto 
regions of low radio continuum emission. 

CO near clump B is probably associated with the IR 
extended source described in Sect. 3.2, as is suggested by 
the striking spatial correlation with the IR ring structure de- 
tected in the far and mid IR (see the bottom panel of Fig. 4). 
Note that the emission in MSX band A is enhanced towards 
the borders of Clump B at (i, b) = (102?21, -0?73). The pres- 
ence of molecular emission in the region indicates that the 
emission in band A most probably originates in polycyclic 
aromatic hydrocarbons (PAHs), which can be found in pho- 
todissociation regions at the interface between ionized and 
molecular gas (e.g. Churchwell et al. 2006). 

We searched for young stellar object (YSO) candidates 
in the IRAS, MSX, and 2MASS point source catalogues 
(Egan et al. 2003, Cutri et al. 2003). Following criteria by 
Junkes et al. (1992), only three IRAS point sources have in- 
frared colors compatible with an YSO classification. Their 
galactic coordinates are listed in Table 2, along with their 
identifications, fluxes in the four IRAS bands, and IR lumi- 
nosities, estimated following Yamaguchi et al. (1999). 

It is worth mentioning that source 2 coincides with 
the extended IR source described above and with clump B. 
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Figure 5. Radio continuum emission at 1420 MHz (contours) superimposed onto the averaged Hi brightness temperature in steps of 3.3 
km s — 1 between -55.5 and —39.0 km s _1 (color scale). The mean velocity of each image is indicated in their upper left corner. Contours 
are 7.4 and 8.0 K. Color scale goes from -15 to 40 K (see text). 



Table 3. H2 column densities and masses. 



Cloud 


W co 


N H 2 


Area 


H2 mass 




K km s — 1 


10 20 cm" 2 


pc 2 


M 


A 


9.9 


10.5±1.4 


14.2 


310 


B 


9.3 


9.9±1.2 


11.1 


200 


C 


6.6 


7.0±0.9 


2.7 


40 


D 


8.9 


9.4±1.2 


12.0 


235 



Source 1 appears projected onto clump A, while source 3 
might be linked to Clump D. The presence of YSO candi- 
dates coincident with molecular clumps suggests that star 
formation is active in certain areas. 



4 DISCUSSION 

4.1 The distance of the ring nebula and the large 
Hi shell 

Since most of the physical parameters of the nebula depend 
on its distance, it is important to estimate its value. Based 
on the association of the WR star to Cep OBI, and adopting 
v = 11.67 mag, and (b—v) = 0.17 mag, van der Hucht (2001) 
places the star at a distance of 2.75 kpc, with an uncertainty 



of about 20%. This distance implies an absolute magnitude 
M v = -2.15 mag for the WR star, somewhat different from 
the adopted mean value for the subclass (in the range -3.9 
to -2.8 mag, see table 24 by van der Hucht [2001]). A quite 
different spectrophotometric distance d* = 4.2 kpc can be 
estimated by adopting a mean value M v = 3.09 mag for 
WN3 stars and the same absorption A v — 1.6 mag (van der 
Hucht 2001). 

As shown in Sect. 3.3, Hi gas related to the ring neb- 
ula around HD 211564 is present in the velocity range from 
—52 to -43 km s _1 , with a mean velocity of ~-47 km s _1 . 
Molecular gas probably related to the inner structure has 
velocities in the range from -49.4 to -46.4 km s~\ Thus, 
a systemic velocity v sya = -47±2 km s _1 for the neutral 
features related to the star seems appropriate. 

The analytical fit to the circular galactic rotation model 
taking into account non-circular motions (Brand & Blitz 
1993) predicts that gas having velocities of -47 km s 
should be placed at a kinematical distance of 3.5 kpc. This 
value is, within errors, compatible with the derived stellar 
distance. As a consequence, we adopt 3.5±1.0 kpc as the 
distance to the WR star and the surrounding structure. At 
this distance, the inner and outer ionized rings have linear 
radii of 8.3±2.4 pc and 15.4±4.4 pc, respectively. 

Hi gas belonging to the large shell has velocities in the 
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Figure 6. Upper panel: Overlay of the CGPS image at 1420 
MHz (contours) and the Hi emission averaged within the velocity 
interval from -47 to —37 km s — 1 . The color scale goes from 
33 to 80 K. No background was removed from the Hi image. 




102.40" 102.20" 102.00" 

Galactic Longitude 



Figure 7. CO(1-0) integrated emission between -56.6 and —46.4 
km s" 1 (contours) superimposed onto the radio continuum emis- 
sion at 1420 MHz (color scale). Contours are 8.2, 14.3, and 20.4 
K km s . 

range -46 to -38 km s -1 , rather similar to the velocity of the 
gas linked to the ring nebula, suggesting a similar distance. 
As a working hypothesis, we adopt a distance of 3.5±1.0 
kpc for the large shell. At this distance, this shell is 31 pc 
in radius. 

4.2 Scenario 

Given the observed morphology, a possible scenario to ex- 
plain the presence of two concentric ionized structures is 
that the inner and outer rings have been formed by the ac- 
tion of the stellar winds from different evolutionary phases 
of the star, as proposed by Garcfa-Segura & Mac Low (1995) 
for NGC 6888. Such a scheme has proved to explain the gas 
distribution in a number of stellar wind bubbles showing 
double structures (e.g. the ring nebulae around WR 16 and 
WR85, Marston 1995). 

However, the presence of neutral atomic gas observed 
between the inner and outer rings, and bordering the inner 
ring, seems to conspire against such an interpretation for 
the case of HD 211564. An alternative explanation is that 
both rings form only one structure, with its main axis nearly 
along the line of sight. This scenario is compatible with the 
detection of gas seen in projection between both rings. 

Here, we propose that the stellar wind bubble associ- 
ated with HD 211564 is evolving in the neutral gas envelope 
of the large shell centered at (l,b) = (102?30, -0?50). Sup- 
port to this interpretation comes from the fact that the outer 
ring near (l,b) = (102? 30, -0?70) appears arched, suggesting 
that it is being pushed from larger galactic longitudes and 
higher galactic latitudes. In this context, the section of the 
outer ring near (102?25, -0?60) would be the result of the 
interaction of both structures. The fact that the ring nebula 
is detected towards the fainter region of the large Hi shell 
and the presence of Hi gas bordering the outer ring at (I, b) 
= (102?00, -0?60) reinforces the interpretation that the in- 
terstellar bubble is evolving in the compressed envelope. In 
this scenario, part of the Hi gas seen inside the inner ring 
and between the outer and inner rings might be linked to the 
large shell. The fact that the large shell and the interstellar 
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bubble have similar velocities makes it difficult to clearly 
identify Hi gas belonging to each individual structure. 

Note that the ionized rings are detected where gas of 
the Hi envelope centered at (102?35, -0?50) is present, i.e. 
where the ambient density is relatively high, in agreement 
with previous findings by Naze et al. (2001) for bubbles in 
the LMC. 



4.3 Physical parameters of the nebula 

We now estimate the main physical parameters of the bubble 
associated with HD 211564 bearing in mind the proposed 
scenario. Thus, we assume that the outer and inner rings 
form only one structure with the radius of the outer ring, 
i.e. 16.3±4.6 pc. 

The parameters of the ionized gas were derived from 
the image at 1420 MHz. Electron densities and Hn masses 
were obtained from the expressions by Mezger & Hender- 
son (1967) for a spherical Hn region of constant electron 
density (rms electron density n e ). The presence of Hell was 
accounted for by multiplying the Hn mass by 1.27. The rms 
electron density and ionized mass are 3.5 cm -3 and 2000 
M©, respectively. Errors of 30% and 60% in the electron 
density and in the ionized mass come from the distance un- 
certainty. 

The distribution of ionized g shown both by the 

optical and radio emissions, suggests filling factors in the 
range / = 0.05-0.12 (estimated assuming a spherical bubble 
of outer radius equal to 16.3 pc, and that 10 to 25% of the 
surface is covered by plasma). By applying these factors, 
electron densities and ionized masses are 10-16 cm -3 and 
450-700 M Q , respectively. 

The number of UV photons necessary to ionize the gas 
in the inner and outer rings, as derived from radio contin- 
uum emission, is log Nl v -c = 48.2, lower than the UV pho- 
ton flux emitted by the WR star (log JV» = 49.2, Crowther 
2007). We can conclude that the WR star can maintain the 
ionization of the rings. The difference between the two val- 
ues is consistent with the fact that a large number of UV 
photons may escape from the bubble through the patchy 
neutral envelope and/or are absorbed by dust grains mixed 
with the ionized gas. 

An important point to discuss in connection with the 
radio continuum emission is the fact that the derived spec- 
tral indices seem to deviate from typical values for thermal 
emission. If the large values of spectral index found for some 
regions (e.g., Fig. [2p) are genuine, this could imply that 
there are regions where the thermal emission is approach- 
ing the optically thick regime. The relative faintness of the 
emission however militates against such an interpretation, 
unless the radio continuum emission originates from regions 
considerably smaller than the observing beam. 

If the emitting regions are in the form of small clumps of 
typical angular size 0, or of thin filaments of angular thick- 
ness 6, beam smearing or dilution will cause the emitted 
regions to have an observed emission measure EM smaller 
than the true emission measure EM by a factor /„ where 
f e « flb/0 2 for clumps and f e ~ 9b /0 for filaments, and fi;, 
and 9b are the beam angular area and beam angular diame- 
ter, respectively. The resolution at 1420 MHz is about 1' so 
that, writing am as the typical angular size in arcmin, we 



simply have f e ~ l/# 2 m an d 1-/0 am for clumps and filaments, 
respectively. 

The brightness temperature due to optically thin brem- 
strahlung or free-free emission Tg is given by Ts = T e th 
where T e and Tg are the electron temperature and free-free 
optical depth, the latter given by (e.g., Mezger & Henderson 
1967, Chaisson 1976): 

T ff = 6.5x 10 lr a^TejT' 1 ' 35 v' 2 ' 1 Eg, (1) 

where a(u, T e ) is the Gaunt factor (about unity for our pur- 
poses) and the emission measure Eg is in pccm -6 . Set- 
ting Tg = 1 in the above equation, and taking a = 1 and 
T e ~ 7000 K, the relation between the turnover frequency 
(in GHz) i^ghz and free-free emission measure is given by 
EM = 1.8 X 10 6 ^G Hz . A turnover at 1 GHz then corre- 
sponds to an emission measure of 1.8 x 10 6 pccm -6 . 

The observed emission measure is given in terms of the 
free-free brightness temperature by EM « 570 Tff which, 
for Tff « 8 K, gives 4500 pc cm" 6 for a = 1 and T e « 7000 K. 
The required beam dilution factor f e is thus of the order 
of 400. To relate this factor to the typical size of clumps 
or thickness of filaments, we note that the linear size of a 
structure in pc is related to its angular size # am in arcmin 
by the relation l pc = 0.29 Sam dk P c, where dk P c is the source 
distance in kpc. For the assumed distance of 3.5 kpc, this 
relation reduces simply to l pc » (9 am . We thus obtain l pc ~ 
1/(400) 1/2 = 0.05 for clumps and l pc » 1/400 = 0.0025 for 
filaments. 

To estimate the neutral atomic mass we took into ac- 
count the Hi gas with velocities between -52 to -43 km s -1 
projected onto the region of the nebula. This velocity inter- 
val corresponds to the range in which Hi gas linked to the 
nebula can be more clearly identified. We assumed that the 
gas is optically thin and included a He abundance of 10%. 
The derived neutral gas mass is 5900 Mq . Uncertainties are 
about 70%. 

The H2 column density and the molecular mass of the 
cloudlets were estimated from the 12 CO data, making use 
of the empirical relation between the integrated emission 
Wco (= fTdv) and N H2 . We adopted N H 2 = (1.06±0.14) 
xWco x 10 20 cm -2 (K km s -1 ) -1 , obtained from 7-ray 
studies of molecular clouds in the Orion region (Digel et al. 
1995). The integrated emission, H2 column density, area, 
and the molecular mass of each cloudlet are listed in Table 
3. 

The small amount of molecular gas in the whole region 
suggests that either little molecular material was present 
in the region when the massive star formed or most of 
the molecular gas was photodissociated and ionized by the 
strong UV stellar flux. The presence of these small molecular 
clouds is compatible with the poor signs of stellar formation 
in the region. 

The ambient density derived by distributing the ionized, 
neutral atomic and molecular mass (5900 M©) over a sphere 
of 16.3 pc in radius is ~14 cm -3 . This value, higher than that 
of the typical interstellar medium, is compatible with the 
suggestion that the interstellar bubble evolved in a neutral 
expanding shell formed by compressed gas. 

To verify whether HD 211564 can provide the energy to 
blow the interstellar bubble, we estimated the mechanical 
energy _E W released by the massive star into the ISM during 
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the dynamical age of the bubble and compared it to the 
kinetic energy Ek of the structure. 

The kinetic energy of the interstellar bubble = 
M b «cx P /2 was derived adopting an expansion velocity of 9 
km s _1 (based on the velocity range of the associated Hi 
gas), and assuming that the total mass related to the bub- 
ble Mb (ionized, neutral atomic, and molecular) amounts to 
5900 Mq. The kinetic energy is 4.8xl0 48 erg. 

The dynamical age of a wind blown bubble can be es- 
timated as t d = 0.56xl0 6 i?/t> OX p yr (McCray 1983), where 
R is the radius of the bubble (pc), w cxp is the expansion 
velocity (km s _1 ), and the coefficient is the deceleration pa- 
rameter. Adopting R = 16.3 pc and v exp = 9 km s~\ id = 
1x10 yr. Large uncertainties are involved in this result. 

The stellar wind mechanical energy E w (— L w t = 
MVyt/2) released by the WR star can be estimated from 
the stellar wind parameters. Since we cannot discard that 
the wind of the O-type star progenitor of the present WR 
star has contributed to the formation of the bubble, we take 
into account the contribution of the star during the WR 
and main sequence phases. The contribution during the WR 
phase can be derived adopting the stellar wind parameters 
listed in Sect. 1 (M = 5x 10~ 6 M Q yr" 1 and V w = 2000 km 
s _1 ), and assuming that this stellar phase lasts 0.5xl0 6 yr, 
compatible with the lifetime of the WR phase of a massive 
star (Meynet & Maeder 2005). E v results to be l.OxlO 50 
erg. To estimate the contribution during the main sequence 
phase, we adopt mean stellar wind parameters for O-type 
stars, i.e. M = 2xl0~ 6 M©yr _1 , V w = 2000 km s" 1 (Prinja 
et al. 1990, Mokiem et al. 2007), and assume that the 
stellar wind acted during at least 3xl0 6 yr. The wind 
mechanical energy released during the main sequence phase 
of the star is then 1.4xl0 5() erg. Consequently, we adopt 
2.4xl0 50 erg as the total energy released through stellar 
winds. 

The energy conversion efficiency e (= E^/E w ) — 0.02 
is compatible with a stellar wind origin for the ring neb- 
ula. Uncertainties in this estimate originate from the input 
stellar wind energy (i.e. the mass loss rates and terminal ve- 
locities), the adopted age of the bubble (which depends on 
the expansion velocity and the expansion law of the bubble) , 
and the kinetic energy (which includes masses and expansion 
velocity) . 

The derived values are similar to estimates for other 
stellar wind bubbles (e.g., Cappa 2006). They are in agree- 
ment with predictions from recent numerical simulations 
from Freyer et al. (2003, 2006). These authors take into ac- 
count the action of the stellar wind and the ionizing flux 
from stars of 35 and 60 M© on the surrounding gas as they 
evolve from the main sequence stage to the presupernova 
phase, and investigate the "missing wind problem". From 
their simulations, they find E^/E w w 0.10-0.04 for stars with 
35 and 60 M© , respectively. Among the additional solutions 
for the "missing wind problem" in superbubbles listed by 
Butt & Bykov (2008), the presence of a blowout, in which 
the hot gas and the energy can spew out from the bubble 
(Cooper et al. 2004) can clearly apply to the present case. 

4.4 The large Hi shell 

Assuming a distance of 3.5 kpc for the large shell, the asso- 
ciated neutral atomic mass amounts to about 10 4 M©. The 



analysis of the CO datacube reveals that a small amount 
of molecular gas may be associated with the shell, indicat- 
ing that the large shell evolved in a region with an ambient 
density of about 3 cm~ 3 . 

The kinetic energy and the dynamical age estimated 
adopting an expansion velocity of 6-9 km s _1 is (4-8) xlO 48 
erg and (1.9-2.6) x 10 6 yr, respectively. 

As regards the origin of the Hi shell, we note that, for 
an energy conversion efficiency of 0.10, an O-type star with 
stellar wind parameters M = 2xlO _6 M©yr _1 and Vw = 
1500 km s _1 might provide the energy to create the ob- 
served structure. However, no O-type star has been detected 
towards the inner part of the large shell. 



5 SUMMARY 

In this paper we have analyzed the distribution of the ion- 
ized and neutral gas, as well as the dust particles, in the 
surroundings of the WR star HD 211564, located at about 
3.5 kpc. 

The radio continuum data shows that the ring nebula 
related to HD 211564 has a clear radio counterpart. Two 
concentric rings of about 9' and 16' in radius are well iden- 
tified at 1420 MHz. The spectral index of the radio con- 
tinuum emission has been estimated in several regions of 
the nebula using the TT-plot method. We have found that 
the spectral index shows significant variations and that it 
slightly deviates from the typical value expected for Hn re- 
gions, suggesting the presence of small and dense clumps 
of ionized gas which are optically thick at 1420 MHz. The 
analysis of the available radio continuum data has enable 
us to estimate some parameters that characterize the Hn 
region. Considering a filling factor in the range 0.05-0.12, 
we obtained ionized masses of 450-700 Mq, and electron 
densities of 10-16 cm~ 3 . 

The analysis of the infrared emission shows that the 
brighter sections of the inner ring are well detected in the far 
infrared, while the outer ring is hardly identified. A strong 
IR source, detected in the four IRAS bands, is observed at 
(l,b) = (102?20, -0?73). This source presents a ring mor- 
phology and is also detected in the MSX Band A. 

An inspection of the Hi images suggests that the cavity 
blown by HD 211564 is evolving in the border of an expand- 
ing shell of about 30' in radius, centred at (I, b) = (102?30, 
-0?50). This large structure is observed in the velocity inter- 
val from -47 to -39 km s _1 . On the other hand, Hi emission 
probably related to the ring nebula is present in the velocity 
range from -52 to -43 km s _ . The fact that both structures 
are partially observed at the same velocities, together with 
their observed morphologies, suggest that both structures 
are interacting. The derived neutral gas mass related to the 
ring nebula is 5900 M©. 

Molecular gas probably linked to the ring nebula is ob- 
served in four small clouds. We have found three YSO candi- 
dates associated with molecular clouds. An energetic anal- 
ysis shows that HD 211564 emits enough ionizing photons to 
keep the gas ionized and heat the dust, and that its stellar 
wind is sufficient to explain the formation of the observed 
features. 
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